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Previous work published at AIAA meeting in Reno 2003, and to be published in AIAA Journal 
of aircraft (FOM=Figure of Merit). 



Motivation? 



Pre-production F/A-18E 

- Exhibited “wing drop” in flight test 

■ “wing drop” is an uncommanded lateral motion 

■ “abrupt wing stall" is an aerodynamic characteristic, and can cause 
wing drop 

- Numerous flight tests resulting in a production fix 

■ Revised flight control laws and porous wing fold fairing 

- A comprehensive program was created to be able to predict these 
phenomenon with wind tunnels and CFD 

< Free-to-roll wind tunnel test method 
- FOM’s for steady and unsteady (non-moving) CFD 

Current work: Progress CFD to calculations of 
damping derivatives and free-to-roll for this flow by 
application to pre-production F/A-I8E 
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For DES, RANS is responsible for predicting boundary layer growth and separation. LES is 
responsible for predicting the geometry dependant turbulent flow features. Grid adaptation done 
using NASA Langley’s RefineMesh program. Adaptation on time average of vorticity 



Previous work (AIAA 03-0594) 


Goal was to use CFD to 
predict the unsteady 
shock oscillations seen in 
the experiments. 

RANS models failed to 
give unsteady results 
Detached-Eddy 
Simulation turbulence 
model 

Hybrid RANS/LES 

- RANS in boundary 
layer 

- LES outside of 
boundary layer 

Solution based grid 
adaption 


Iso of vorticity, colored by pressure, a=9" 
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DES results are time-averaged coefficients. Left axis removed to protect proprietary data. 


a b a I t 


Previous work (AIAA 03-0594) 


10/10/5 flap set with no 
tails 

SST predicted early lift 
curve break 

DES showed ail improved 
lift curve break (but on a 
grid finer than the current 
grid) 

Motivates inclusion of 
DES in the current project, 
along with RANS 
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These DES projects represent a cross section of those done over the past few years using Cobalt. 

Delta wing vortex breakdown on a delta wing and the F-18C done by Major Scott Morton of the 
USAF Academy (Scott.morton@usafa.af.mil). 

2-D forebody geometry by Kyle Squires (squires@asu.edu). 

Prescribed spin of the F-15E by James Forsythe. 



438 


Prisms created using “Blacksmith” to recombine the tets in the boundary layer into prisms. 
Blacksmith is a Cobalt grid utility. 




Grid mirrored about symmetry 
plane 

Grid provided by Paresh Parikli 
6/8/4 flap set 

8.4x1 0 6 cells for both sides of 
aircraft 

Adaption performed 
on a 9° time-averaged 
DES solution under 
previous work 

Prisms in boundary layer 
Average y + < 0.7 


Grid 
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The following are non-moving cases - but can be unsteady (for DES) 
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CPU hours based on a Compaq ES45. Timestep for DES non-dimensionalized by chord and 
freestream velocity. 


) 

Solution Procedure 

Menter's SST RANS model 

- Convergence monitored by observing forces and moments. 
Rolling moment was generally the most sensitive and last 
to converge. 

* 4000 iterations 

* 1 Newton sub-iteration 
. CEL of lxl 0* 

* 2000 epu hours per run. M w 

Spalart-Allmaras based DES model 

- Unsteady flow simulation 

• 16000 iterations 

♦ 3 Newton sub-iterations 

\ 

• 8x the cost of the steady RANS simulations 
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Model was set to a given pitch angle (theta), then rolled about the longitudinal axis (phi). This 
resulted in a decrease of alpha, and an increase in beta as phi increased. The CFD was 
performed at the given alphas and betas, which were corrected in the wind tunnel data for wall 
effects. 


Test Matrix 


Conditions chosen to 
match NASA Langley 
wind tunnel test 565 

Mach=0.9 

Re c =3.9xl0 6 

Flow through engines 

Sting not included in 
grid 


ALPHA 

BETA 

PHI 

THETA 

7.11 

-0 14 

-0 13 

7 01 

7 03 

0.86 

9 87 

6 97 

6 32 

2 73 

29 75 

6 74 

4 03 

4.89 

59.78 

6.14 

240 

ALPHA 

BETA 

PHI 

THETA 

B 50 

0.35 

3 84 

8 40 

8 12 

1 04 

9 82 

3 08 

7 33 

3 42 

29 63 

8 18 

242 

ALPHA 

BETA 

PHI 

THETA 

9 04 

1.17 

9.69 

9 00 

8 20 

3 74 

29 56 

8 87 

244 

ALPHA 

BETA 

PHI 

THETA 

10 17 

-0 12 

-0.02 

10.06 

10.05 

1 38 

9.79 

10.03 

9 06 

4.23' 

29 49 

9 86 

5 69 

7.54 

59 53 

9.26 
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CNFBA 



Left axis not labeled to protect proprietary data 



Narmal force for near-zero sideslip 
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CPMBA 



Pitching moment for near-zero sideslip 
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CRMBA 


Note reversal of rolling moment for phi=30 using SST. 
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CYMBA 


Yawing moment well predicted - as with all cases. 



Run 247 (0*7°) 
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CSFBA 


Side force well predicted - as with all cases. 
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Shock retreating off trailing edge of leading edge flap. 
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DES isosurface looks like separation is at trailing of leading edge flap. But it moves back from 
there unsteadily. This leads to the blue low pressure in the separation bubble (since it is not 
always separated). 


I so surface of u=0, surface colored by pressure 


Run 247 (0«7°), (j)«0° 

SST DES 
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Isosurface of u=0, surface colored by pressure 


Run 247 ( 0«7 ), (|)«10 

SST DES 
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The separation moving forward on the right wing is the cause for the roll moment reversal. 
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At this high phi, the alpha is reduced so much that the flow remains attached until the trailing 
edge of the wing. 
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GRMBA 


Note asymmetries in wind tunnel data. Decrease in lateral stability derivative picked up with 
DES. 



Run 240 (0«8.5°) 
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CYME A 


Good agreement for yawing moment, as with all cases - this is likely due to the attached flow at 
the tail, which is easily predicted. 



Run 240 (0*8.5°) 
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CSFBA 


Good agreement for side force, as with all cases - this is likely due to the attached flow at the 
tail, which is easily predicted. 


Run 240 (0«8.5°) 


Side Force 
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Separation is making it onto the leading edge of the leading edge flap. 
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Isosurfacc of u-0, surface colored by pressure 

Run 240 (0*8.5°), <j>»4° 

SST DES 
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Isosurfacc of u=0, surface colored by pressure 


% Run 240 (0*8.5°), (b*10° 


DES 


SST 
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Isosurface of u=0, surface colored by pressure 

Run 240 (0*8.5°), c|>*30 o 
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CRM6A 


Large asymmetries in wind tunnel data. Around this angle there was difficulty in testing, since 
model dynamics became significant. 
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CYM6A 


Good agreement for yawing moment, as with all cases. 
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CSFBA 


Good agreement for side force, as with all cases. 



Run 242 (0«9°) 
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CRMS A 


Rolling moment offset predicted by DES - is the sample size large enough? 
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Looks like enough samples have been taken to well define rolling moment. However more 
might change the time-averaged rolling moment some. 
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467 
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Unsteadiness now is due to separation moving from leading to trailing edge of the leading edge 
flap. 



JpBS, Run 244 (0*1 0°), 4>«0° 

Isosurfacc of vortlcity, colored by pressure 
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Isosurface of u=0, surface colored by pressure 


Run 244 (0*1 0 ), (|)«0 

SST DES 
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I sd surface of u=0, surface colored by pressure 

Run 244 (0*1 0°), 4>«1 0° 

SSF DES 
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ALE = Arbitrary Eulerian/Lagrangian. 





Solution Procedure 


Time-accurate with ALE formulation for grid motion 

- 5 Newton sub-iteration (for accurate grid motion) 

- At*=0.02 (ran several timesteps to demonstrate timestep 
convergence) 

Prescribed sinusoidal oscillation around longitudinal 
axis 

- tan _1 (f ,s )=l 0 ; f*=0.0174 

- 2,600 iterations per cycle 

- 4,000 epu-hours per cycle 

+/-5° oscillation v 

Menter’s SST RANS model 


= tan 
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Linear and well behaved. Stable roll damping. 
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Separation at trailing edge - flow well behaved. 



All 



Large rolling moment offset. Several cycles run with varied timestep, but offset remained. 


oiling moment vs. Roll rate 



cty/cT b/2U 
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Offset due to differences in separation location. Hysterisis? 
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Slightly chaotic behavior, but linear and stable roll damping. 
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Positive roll damping. Note lowered slope - due to lower lift curve slope once shock moves 
forward on the wing. 



0.02 

d4>/dt*b/2U ; 


0 = 5 “ 

6 = 6 “ 

6 = 7 ° 

6 = 8 “ 
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oiling moment vs. Roll rate 
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oiling moment vs. Roll rate 
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Roiling Moment vs. Roll angle 
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Study still underway - not enough samples. Looking at dependence of roll damping on roll rate. 



Roll rate study, 0=7° 


atan(f)=Q.5“ 

alan(f)-1 n 

atan(f)=2" 


o 










m 


-J L L- 


" s '~i 


- 0.04 - 0.02 0 0.02 

d4>/dt*b/2U 


0.04 
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Conclusions 

RANS and DES applied to predict static stability 
derivatives in roll in AWS regime 

- DES showed better lift and moment predictions 

- Yawing moments and side force well predicted by both methods 

- Rolling moment more sensitive (both for CFD and wind tunnel) 

Prescribed rolls used to look at roll damping (RANS only) 

- All cases were stable in roll, but in AWS regime had more 
chaotic behavior. For one angle there was a significant rolling 
moment offest 

- Comparison to experiments still ongoing 

• Continuing work 

- DES of prescribed rolls 

- More iterations on varying roll rate 
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